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Using a microscopic transport model we investigate the evolution of conical structures originating 
from the supersonic projectile moving through the hot matter of ultrarelativistic particles. Using 
different scenarios for the interaction between projectile and matter, and different transport proper- 
ties of the matter, we study the formation and structure of Mach cones. Especially, a dependence of 
the Mach cone angle on the details and rate of the energy deposition from projectile to the matter is 
investigated. Furthermore, the two-particle correlations extracted from the numerical calculations 
are compared to an analytical approximation. We find that the propagation of a high energetic 
particle through the matter does not lead to the appearance of a double peak structure as observed 
in the ultrarelativistic heavy-ion collision experiments. The reason is the strongly forward-peaked 
energy and momentum deposition in the head shock region. In addition, by adjusting the cross 
section we investigate the influence of the viscosity to the structure of Mach cones. A clear and 
unavoidable smearing of the profile depending on a finite ratio of shear viscosity to entropy density 
is clearly visible. 



I. INTRODUCTION 

Results from the relativistic Heavy Ion Collider 
(RHIC) p] and recently from the Large Hadron Collider 
(LHC) [2 indicate the formation of a new state of matter, 
the quark-gluon plasma (QGP). The large value of the 
measured elliptic flow coefficient i>2 indicates the nearly 
perfect fluid behavior of the QGP 3 . This is confirmed 
by recent calculations of viscous hydrodynamics [5j and 
microscopic transport calculations [5 with a shear viscos- 
ity over entropy density ratio r//s — 0.1 — 0.2, which is 
close to the conjectured lower bound Tj/s = 1/Att from a 
correspondence between conformal field theory and string 
theory in an Anti-de-Sitter space [B]. 

Highly energetic partons propagating through the hot 
and dense QGP rapidly lose their energy and momen- 
tum as the energy is deposited in the medium. This 
phenomenon is known as jet-quenching [5], whereas 
its exact mechanism is still to be fully understood. Fur- 
thermore, recent measurements of two- and three-particle 
correlations in heavy-ion collisions (HIC) show a com- 
plete suppression of the away-side jet, whereas for lower 
Pt a double peak structure is observed in the two-particle 
correlation function [9 . For a while one possible and 
promising origin of these structures was assumed to be 
the interaction of fast partons with the soft matter which 
generates collective motion of the medium in form of 
Mach cones. [10l[TT]. In contrast, recent studies of tri- 
angular flow from initial fluctuations |12j show a more 
satisfactory explanation for the appearance of the dou- 
ble peak structure. 

The recent idea is that both Mach cones and triangular 
flow from initial fluctuations exist in heavy-ion collisions, 
but it is difficult to separate their effects and the two- 
particle correlations do not seem to be a good observable 
for this purpose. In the present study we claim that even 



if there are no effects from initial stage fluctuations, a 
double-peak structure in two-particle correlations cannot 
be expected from the energy-momentum deposition by a 
jet into the medium. This will be not only due to the 
viscous effects, but also due to the details of the "Mach 
cone" -like structure, which was found in a similar form 
in ideal fluid [13] and AdS/CFT[14 studies. 

For this purpose we investigate the propagation and 
formation of Mach cones in the microscopic transport 
model BAMPS (Boltzmann Approach of MultiParton 
Scatterings) |15j in the limit of vanishing mass and very 
small shear viscosity over entropy density ratio i}/s of 
the matter. Two different scenarios for the jet are used 
and the dependence of the Mach cone angle on the de- 
tails of energy deposition is discussed. A simple ana- 
lytic relation for the expected particle distribution in the 
Mach cone wings is derived and is compared it to nu- 
merical results extracted from BAMPS. In addition, by 
adjusting ?y/s, the influence of the viscosity on the pro- 
file of the Mach cone and the corresponding two-particle 
correlation is explored for the first time. In this work 
the units are % = c = k = 1. The metric tensor is 
cT =diag(l, -1,-1,-1). 



II. SHOCK WAVES AND MACH CONES 

Shock waves are phenomena which have their origin in 
the collective motion of matter. In the limit of a perfect 
fluid with no viscosity, a signal caused by a weak pertur- 
bation propagates with the speed of sound c s = \J dp/de, 
which depends on the equation of state (EOS) of the 
medium. Here, p is the equilibrium pressure and e is the 
energy density in the local rest frame (LRF). A larger 
perturbation results in a shock wave propagating faster 
than the speed of sound. In a simplified one-dimensional 
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setup shock waves have already been studied for the per- 
fect fluid limit [16 . Furthermore, the viscous solutions 
have been investigated in Refs. [TTJ [JJ5], demonstrating 
that the shock profile is smeared out when viscosity is 
large. It was also found that a clear observation of the 
shock within the short time available in heavy-ion colli- 
sions requires a small viscosity. The information taken 
from these studies can be transferred to the investigation 
of conical shock structures like Mach cones, which is the 
main subject of this Letter. 

First, we consider a weak perturbation moving with 
the speed of light, i.e. u SO urce = 1, through the medium 
of a perfect fluid. For simplification we use a massless 
relativistic gas with e = 3p and c s = l/v3- The pertur- 
bation generates waves propagating through the medium 
with the speed of sound c s . In this case the propagating 
modes are called sound waves. If the perturbation moves 
faster than the speed of sound the created sound waves 
accumulate on a cone [T5]. In the following we refer to 
the surface of this cone as the shock front. 

The resulting emission angle of this shock front rel- 
ative to the direction of the projectile is given by the 
weak perturbation Mach angle a w = arccos(c s /u source ) = 
arccos(l/-\/3) = 54,73°. It is important to know that in 
nature perturbations are not sufficiently small. In this 
case, shock waves instead of sound waves are generated 
and due to different propagation velocities of these waves, 
we expect a change of the Mach angle [5D]. We can gen- 
eralize the introduced Mach angle to the case of stronger 
perturbations: 



arCCOs(t> s hock/^sourcc) • 



(1) 



We require here Source > v s hocki where u s hock is the veloc- 
ity of the shock front propagating through the medium. 
The velocity of the shock front depends on the pressure 
(energy density) on the conep cone (e con e) and the medium 
itself p med (e med ) [IB]: 
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Eq. ([2]) has the following limits: If p CO nc ^> p m cd we ob- 
tain w s hock = 1- If Pconc ~ Pmcd, i-e. the perturbation 
is very weak, we get the expected limit of the speed of 
sound Wshock ~ c s . In the latter case Eq. (JlJ becomes 
expected. The collective velocity of matter in the 
shock wave (Mach cone wing) , which is different from the 
signal propagation velocity |2|, can be calculated via 

(Pconc Pvacd ) (^conc ^-mcd) / \ 
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In the case of a very weak perturbation the collective ve- 
locity of matter vanishes, v co ii ~ 0; whereas for stronger 
perturbations i? co u can increase up to the speed of light. 

In Sec. [TV] we discuss the numerical results from 
BAMPS and expect a clear dependence of the observed 
Mach angle on the strength of the perturbation according 
to Eq. , but due to non-linear effects Eq. is merely 
a good approximation. 



III. PARTICLE MOMENTUM DISTRIBUTION 
IN THE SHOCK FRONT 



In order to understand the origin of the double peak 
structure induced by "Mach cone" -like structures, which 
will be discussed in Sec. |IV| we derive a simple model of 
particle emission from the shock front of a Mach cone in a 
2-dimensional xy-plane. We assume two sources model- 
ing the two wings of a Mach cone with a constant temper- 
ature T and collective four- velocity = 7(1, v), where 
7 = — v 2 is the Lorentz gamma factor. Each source 

consists of massless particles according to the thermal 
distribution f(x,p) — exp(— u^p^/T), where p^ = (E,p) 
is the particle four-momentum. Choosing the x-axis to 
be the symmetry axis of the cone, which is simultane- 
ously the propagation direction of the jet, we can write 
u± = 7(1, v cos a, ±1; sin a, 0). The ± corresponds to 
each wing of the cone. We identify v — v co n with Eq. ^ 
as the collective velocity of the matter in the shock wave 
and a is the Mach angle defined in Eq. ([I]). Using the 
same coordinate system we write for the four-momentum 
vector p^ = p(l, cos cf> sin 9, sin (f> sin 9, cos 9). 4> is the az- 
imuthal angle in the xy plane and 9 is the polar angle 
with the z-axis. 

The distribution function is defined as 
dN(2irf/{dVd 3 p) = f(x,p), where dVd 3 p/(2Tr) 3 
is the phase space volume element. We are interested 
in the particle distribution dN/ (Nd(f>) which can be 
calculated as an integral over the thermal distribution in 
a certain volume V on the Mach cone surface. We use 
d 3 p = p 2 dp d(j> rf(cos 9) and write 
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We obtain N = 8-K-fT 3 V by integrating dN/(dVd 3 p) = 
f(x,p) over the entire phase space volume. After the 
integration of Q we obtain 
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where A = tt/2 + arctan(6j/y / l — b 2 ), b\ — vcos(a — <$>) 
and 62 = ucos(a + 4>). 

The most important result taken from Eq. |5]) is the 
non-existent double peak structure for small w C oib which 
is against all expectations resulting from the naive pic- 
ture of a Mach cone. We will discuss this point in more 
details in Sec. IV by comparing Eq. ([5| to the numerical 
results. Furthermore, the particle distribution is equiva- 
lent to the two-particle correlation, since the angle <f> is 
always correlated to the direction of the source, which 
serves as a "trigger" particle. 
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IV. TRANSITION FROM IDEAL TO VISCOUS 
MACH CONES IN BAMPS 

In the following we study the evolution of "Mach 
cone"-like structures with different scenarios of the jet- 
medium interaction by using the parton cascade BAMPS 
[15j - a microscopic transport model which solves the 
Boltzmann equation p ,1 d^f{x,p) = C[f(x,p)] for on- 
shell particles based on stochastic interpretation of tran- 
sition rates. As was demonstrated in previous works, 
BAMPS is able to explore a large variety of hydrody- 
namic phenomena and provides a reliable benchmark 
for hydrodynamic models [17l [21]. The advantage of 
BAMPS is its ability to handle arbitrary large gradi- 
ents for any choice of viscosity. Thus it is possible to 
investigate the complete transition from ideal to viscous 
behavior. 

In this study we focus on investigation of Mach cone 
evolution in absence of any other effects - i.e. we ne- 
glect such effects as initial fluctuations or expansion, 
which are however relevant in heavy-ion collisions. For 
this purpose, the space-time evolution of particles is per- 
formed in a static box. We initialize a static uniform 
medium of massless Boltzmann particles with T me( j = 400 
MeV, which corresponds to a LRF energy density e mc( j = 
16.28 GeV/fm 3 . For simplification, we consider only bi- 
nary collisions with an isotropic cross section, i.e. a cross 
section with an isotropic distribution of the collision an- 
gle. Furthermore, we keep the mean free path A m f p of the 
medium particles constant in all spatial cells by adjust- 
ing the cross section according to a — l/(nA m f p ), where 
n is the particle density. The related shear viscosity for 
isotropic binary collisions is given by r\ = 0.4eA m f p [2"2"] . 
Collisions of particles against box boundaries in x and 
y direction are realized as elastic collisions off a wall; in 
z-dircction we use periodic boundary conditions. This 
reduces the problem to two dimensions and therefore de- 
creases the numerical expenses. 

We introduce two different sources to investigate the 
evolution of "Mach cone" -like structures. In the so called 
pure energy deposition scenario (PED) |13j the source 
propagates and emits particles according to the thermal 
distribution f(x,p) — exp(—E/T), so that the energy de- 
position is isotropic in the LRF of the source. In this sce- 
nario on average only energy is deposited to the medium, 
but no net-momentum. In the second scenario, referred 
to as JET, a highly energetic massless particle (jet) has 
only momentum in x-direction, i.e. p x = Ej Ct . The jet 
propagates and deposits energy to the medium due to 
collisions with particles. After each collision, the mo- 
mentum of the jet is reset to its initial value. The jet- 
medium cross section is adjusted in such a way that we 
obtain a specific energy deposition rate. Using this sce- 
nario a constant energy and momentum deposition rate 
is achieved. For both scenarios the sources are initialized 
at t = fm/c at the position x = —0.1 fm and propa- 
gate in x-direction with i> source = 1, i.e. with the speed 
of light. We note the JET scenario is a simplified model 



of a jet in heavy-ion physics, whereas the PED scenario 
vaguely resembles the hot spots studied in [23], but in 
the form implemented here there is no correspondence to 
heavy-ion collisions. We expect clear differences between 
these two scenarios concerning the evolution of the entire 
system, but also concerning the final distribution of the 
particles. 



A. Effect of energy deposition rate 

In Fig. [T] we show the results for the PED in the upper 
panel and JET in the lower panel using three different 
energy deposition rates into the medium, dE/dx = 1, 10 
and 200 GeV/fm, in the nearly ideal limit, i.e. r\j s k, 

0. 005. We note that in general the maximum (minimum) 
energy density in the simulations is larger (smaller) than 
the maximum (minimum) of the energy density scales in 
the figures. Also the plotted arrow length of the velocity 
profile is scaled. Both modifications are done to enhance 
the readability of the figures. 

In both scenarios, PED and JET, we observe a coni- 
cal structure, but with obvious differences. In the PED 
case with the isotropic energy deposition, a circle of per- 
turbations propagating in backward direction is visible. 
This is missing in the JET scenario because of the strong 
momentum deposition in x-direction. Another difference 
is that in the JET scenario a clearly visible head shock, 

1. e. a shock wave in the front of the jet perpendicular to 
the direction of the jet, appears. This in turn is missing 
in the PED scenario. Furthermore, there is a clear dif- 
ference in the behavior of the matter behind the Mach 
cones. In the JET case, the projectile induces a diffu- 
sion wake, where the matter is flowing in the direction 
of the projectile. Whereas in the PED scenario an oppo- 
site behavior is observed, i.e. there is an (anti-) diffusion 
wake where the matter behind the cone is flowing in the 
backward direction. These observations are in qualitative 
agreement with the results from ideal hydrodynamics and 
transport calculations [TBI 12"4"] . 

Additionally, every scenario is compared to the ideal 
Mach cone with a w for a very weak perturbation shown 
m Fig. [1] Both scenarios provide evidence that the en- 
ergy deposition rate of the source influences the Mach 
angle a of the wings according to Eq. [T] In both cases 
the shock front is curved because near the projectile the 
disturbance of the media is strongest and the shock front 
moves faster than the speed of sound. Farther away from 
the projectile a part of the energy of the shock front has 
already dissipated into the medium and as a result the 
perturbation gets weaker and approaches a weak pertur- 
bation propagating with the speed of sound. 

In the JET scenario the energy of the jet £j t is 20, 
200 and 20000 GeV (starting from the left in Fig. [I}. 
For our calculations in the nearly ideal limit the energy 
of the jet does not play any significant role. The only 
parameter which matters is the average energy deposition 
rate. We will mention in Sec. HVBI how the value of the 
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FIG. 1. (Color online) Shape of a Mach cone in the nearly ideal limit (77/s = 0.005) shown for different jet scenarios and different 
energy deposition rates into the medium, dE/dx = 1, 10 and 200 GeV/fm. The upper panel shows the pure energy deposition 
scenario (PED); the lower panel shows the propagation of a highly energetic jet (JET) depositing energy and momentum in 
^-direction. Depicted are the LRF energy density within a specific range; as an overlay we show the velocity profile with a 
scaled arrow length. The results are a snapshot of the evolution at t — 2.5 fm/c. In addition we show the analytical solution 
for the ideal Mach cone in the very weak perturbation case with the emission angle a w . 



jet energy E- ]et and finite viscosity in the medium changes 
the pattern of the Mach Cone. 

We now want to address the question whether the 
Mach cone structures observed in Fig. [T] can be regarded 
as the source of a double peak structure in two-particle 
correlations. For this purpose we extract the particle 
distribution dN/(Nd~4>) from BAMPS calculations. In 
Fig. [2] (a) we show the results for the energy deposition 
rate dE/dx — 10 GeV/fm together with the analytical 
calculation using Eq. ([5]). To extract only the contri- 
bution from the wings and to exclude of all other re- 
gions such as (anti-) diffusion wake and back region (es- 
pecially in the PED scenario) , a lower energy density cut 
at 20 GeV/fm 3 is applied. Particles in cells with energy 



density lower than this value are not considered in the ex- 
tracted particle distribution (we note that particles from 
the medium in rest automatically do not contribute to 
the final profile). For the analytical solution taken from 
Eq. 45) we use e conc = 22.15 GeV/fm 3 and v co n — 0.137 
(e CO ne represents the average energy density on the Mach 
Cone wings extracted from the associated numerical cal- 
culations). In both scenarios, PED and JET, as well as 
in the analytical calculation we observe only a peak in 
the direction of the source, but no double peak struc- 
ture. This finding is against all expectations from the 
naive picture of a Mach cone. 

However, with a sufficiently higher energy deposition 
rate the final picture changes significantly. In Fig. [2] (b) 
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FIG. 2. (Color online) Two-particle correlations dN/(Nd(f>) 
extracted from calculations shown in Fig. [l] The results are 
extracted from calculations with dE/dx = 10 GeV/fm (a) 
and dE/dx — 200 GeV/fm (b). Analytic solutions extracted 
from Eq. ([5| are shown for e cono = 22.15 GeV/fm 3 (a) and 
e cone = 62.55 GeV/fm 3 (b). 



the results from BAMPS calculations with dE/dx = 200 
GeV/fm are shown. The lower energy density cut is 
increased to 50 GeV/fm 3 because of the much higher 
energy deposition rate. For the analytic calculation 
econe = 62.55 GeV/fm 3 with v co \\ = 0.537 is selected. 
In the PED scenario, as well as in the analytic model, 
the double peak structure finally appears as long as the 
energy deposition rate and consequently v co \\ are suffi- 
ciently large. However, in the JET scenario only a peak 
in the direction of the jet is visible. 

We want to mention that in the PED scenario using 
special momentum cuts (not considered in this work), 
i.e. restricting the momentum integration in Eq. ^ to a 
certain interval one can always obtain a double peak in 
the distribution. However, in the JET scenario a double 
peak never appears, regardless of momentum cuts. 

There are two main contributions to the structure of 
the two-particle correlation, one from the wings of the 
Mach cone and one from the head shock region. The mat- 
ter in the wings is moving in the direction perpendicular 
to the surface with some collective velocity v co \\. The 
larger the collective velocity, the more strongly peaked 
are the local particle distribution functions into this di- 
rection. From our simple analytic model it is clear that 
the mere existence of the wings is not enough to have 
clearly visible peaks in the correlation, but the local ve- 
locity of the matter has to be sufficiently large. This is 
also confirmed by the full simulations: If the energy de- 
position rate is sufficiently large in the PED scenario, the 
double peaks appear. 

In principle, the same reasoning also works for the JET 
scenario. However, in this case there is also a strong con- 
tribution from the head shock region, where the matter 
is moving with large collective velocity. This collective 
motion is in the direction of the projectile and results 
in a particle distribution function that is peaked in the 
same direction. Although a double peak due to the Mach 
cone wings still exists, the contribution of the head shock 
clearly dominates and overshadows the contribution from 
the wing regions (with spatial cuts to remove the head 
shock the double peak appears again [TT]). Thus, no 



double peaks appear in the JET scenario. 



B. Effects of viscosity 

In Fig. [3] we show the Mach Cone structure for both 
PED scenario (upper panel) and JET scenario (lower 
panel) with rj/s = 0.005, 0.05, 0.2 and 0.5 from left to 
right, respectively. The energy deposition rate is fixed to 
dE/dx = 200 GeV/fm. In addition, E }et = 20000 GeV 
is used in the JET scenario. The chosen rj/s values are 
intended to cover the nearly-ideal limit (0.005), the es- 
timated QGP shear viscosity over entropy density ratio 
in heavy-ion collisions (0.05, 0.2) [3] [5] and highly viscous 
limit where dissipative hydro calculations are not reliable 
anymore (0.5) [T7] . 

First, we note that if we observe the system at fixed 
time, then in both scenarios the Mach cone structure 
smears out and eventually disappears almost completely 
as the viscosity increases. This is true for shock fronts 
as well as for the (anti-) diffusion wake. The difference 
between the PED and the JET case is that as rj/s in- 
creases, in the PED scenario the resulting "Mach cone" 
solution covers approximately the same spatial region re- 
gardless of a value of rj/s, while in the JET case the 
structure is concentrated more and more near the pro- 
jectile as the viscosity increases. The reason for this is 
that in the PED scenario the momentum from the pro- 
jectile is isotropically deposited into the medium, while 
in the JET scenario the initial momentum dissipation is 
strongly peaked into the direction of the projectile (the 
effect in the JET scenario becomes even stronger with in- 
creasing energy of the jet Ej ct , since scattered the parti- 
cles are stronger forward-peaked) . With a large viscosity 
the re-scattering of the emitted particles from the source 
is very rare. Thus, the larger the viscosity the more the 
resulting solution reflects the details of the projectile- 
matter interaction. 

We note that in both scenarios the projectiles are 
point-like and initially the matter is homogeneously dis- 
tributed. Therefore, as long as we keep the energy dissi- 
pation rate constant, the only length scale that controls 
the solution is the mean free path A m f p oc rj. Thus, we ex- 
pect a similar scaling behavior as in the one-dimensional 
Riemann problem |17) . For example, the energy density 
profiles for two different shear viscosities rj and rj are 
related by 

/. , i\ f t -to x-x Q y-y , 

e(t-t ,x-x ,y-y ,ri) =e{—^, - , - ,77), 

(6) 

where the scaling factor C = rf /rj, and Xq and yo are the 
coordinates of the projectile at the time to- 

Using this scaling behavior, we can also read Fig. [3] as 
a time-evolution of the solution, with a larger viscosity 
corresponding an earlier time. For example, the solu- 
tions with rj/s — 0.5 in the right-most panel of Fig.[3]will 
evolve to the ones with rj/s — 0.05 at time t = 25 fm/c 
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t = 2.5 fm/c; dE/dx = 200 GeV/fm 
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FIG. 3. (Color online) Transition from ideal to viscous Mach cones. Shape of a Mach cone shown for different jet scenarios and 
different viscosity over entropy density ratios, rj/s — 0.005, 0.05, 0.2 and 0.5. The energy deposition is dE/dx = 200 GeV/fm. 
The upper panel shows the pure energy deposition scenario (PED); the lower panel shows the propagation of a highly energetic 
jet (JET) depositing energy and momentum in :r-direction. Depicted are the LRF energy density within a specific range; as an 
overlay we show the velocity profile with a scaled arrow length. The results are a snapshot of the evolution at t — 2.5 fm/c. In 
addition we show the analytic solution for the ideal Mach cone in the very weak perturbation case with the emission angle a w . 



(with the appropriate scaling of the x— and y— axis). Al- 
though, from Fig. [3] the Mach angle apparently changes 
with the viscosity, this is a transient effect related to a 
finite formation time of the Mach cone with non-zero vis- 
cosity. The viscosity affects the width and formation time 
of the shock front, but not its speed of propagation, i.e. 
the relation (|2| still holds for non-zero viscosity. Asymp- 
totically, the Mach cone angle will be the same regardless 
of the value of rj/s. 
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FIG. 4. (Color online) Two-particle correlations dN/(Nd(f>) 
for different viscosities extracted from calculations shown in 
Fig. [3] The results are shown in the for the JET (a) and PED 
(b) scenario for dE/dx = 200 GeV/fm. 

In Fig. [4] we show the two-particle correlations for the 
solutions from Fig. [3] The procedure is similar to the 



one discussed for Fig. [2j The lower energy density cut is 
chosen to be 50 GeV/fm 3 . For the JET scenario (a), the 
peak in direction of the jet becomes sharper with larger 
viscosity and no other appreciable effect originating from 
viscosity is visible. In contrast, for the PED scenario (b) 
the viscosity destroys the double peak structure. If the 
viscosity is very large, only a peak in direction of the jet 
is visible. As above, these results can also be read as 
a time-evolution of the solution. Fig. [4] shows how the 
angular distribution of the emitted particles widens with 
time, or equivalently with increasing rj/s. 



V. CONCLUSIONS 

In summary, we have investigated the structure of rel- 
ativistic Mach cones by using a microscopic transport 
model. The simulations were realized by using two dif- 
ferent types of sources propagating through the matter, 
the PED and JET scenarios. The effect of the strength 
of the projectile-matter interaction was studied by vary- 
ing the energy dissipation rate from the projectile to the 
matter. Furthermore, the effect of the viscosity of the 
matter was investigated by adjusting the shear viscosity 
over entropy density ratio rj/s from 0.005 to 0.5. 

We observed the conical structures form for both types 
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of sources in the nearly perfect fluid limit, similar to ob- 
servations in [13], with the Mach Cone angle depending 
on the energy dissipation rate. We also demonstrated 
that the non- vanishing viscosity tends to destroy the clear 
conical structure. By using a scaling of the solutions, we 
argued that increasing the viscosity has the same effect 
as looking at the solution at an earlier time. The larger 
the viscosity or, equivalently, the less time the Mach cone 
has to develop, the more the structure of the solution de- 
pends on the details of the projectile-matter coupling. 

Although Mach cone-like structures are observed in 
BAMPS calculations for different energy and momentum 
deposition scenarios they are not necessarily associated 
with double peak structures in the azimuthal particle dis- 
tributions in dN/(Nd4>). We found that only the PED 
scenario together with a rather high energy deposition 
rate lead to a double peak structure, which otherwise 
cannot be observed because of the strong diffusion wake 
and head shock. However, the PED scenario has no cor- 
respondence in heavy-ion physics. On the other hand, 
the JET scenario is a simplified model but neverthe- 



less demonstrates that a double peak structure cannot 
be produced by jets with energy and momentum depo- 
sition. We expect that our conclusions will still be valid 
for realistic jets and energy loss scenarios [25] . In addi- 
tion, a clear Mach cone structure, which is necessary but 
not in itself sufficient to produce a double peak structure 
in two-particle correlations, will hardly develop in a sys- 
tem of the size and finite viscosity relevant for HIC. We 
thus conclude that the double peak structure is not the 
appropriate observable for the signal of Mach cones in 
heavy-ion collision experiments. 
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